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It  is  well  known  that  morphological  effects  play  a  very  important  role 
in  surface  spectroscopic  studies  [1],  and  in  particular,  decay  rates  for 
molecules  located  near  a  sphere  (which  serves  as  a  prototype  surface 
structure)  have  been  reported  extensively  in  the  literature  [2-6]. 

Previously,  the  behavior  of  the  decay  rates  as  functions  of  the  molecule- 
surface  distance  and  molecular  orientation  has  been  studied.  Tn  addition, 
it  has  been  observed  that  the  simple  nonretarded  image  theory  can  be  very 
inaccurate  for  large  molecule -surface  distances  or  highly-conducting  spheres 
[5].  In  this  work,  we  shall  study  further  the  decay  rate  spectrum  as  a 
function  of  emission  frequency  within  the  context  of  the  dynamic  energy 
transfer  theory  (ET)  [3-6]  and  compare  it  with  the  results  obtained  from  the 
static  image  theory  (IT)  [2].  We  shall  see  that  the  morphology-dependent 
peak  obtained  previously  by  Ruppin  [3]  can  be  interpreted  as  the  radiative 
energy  transfer  resonance  between  the  molecule  and  the  sphere,  and  that  such 
peak  can  be  described  only  in  the  context  of  the  dynamic  theory . 

2 .  Brief  Summary  of  Previous  Results 

Let  us  consider  the  simplest  case  of  a  radial  dipole  located  at 
(d,0,0)  above  a  sphere  of  radius  a.  By  assuming  a  quantum  yield  of  unity, 
the  total  (T)  and  radiative  (R)  normalized  decay  rates  have  been  obtained 
previously  within  ET  in  the  form  [6] 
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where  j  and  h^^  are  the  spherical  Bessel  and  Henkel  functions,  and  k  -  w/c 

is  the  wave  number  of  the  molecular  emission.  B  is  a  kind  of  "Mie 

n 

coefficient"  defined  by 


jn(ka)  ^(/7ka)  -  e  jn(77ka)  V>^(ka) 
n  <  jn(77ka)  f;(ka)  -  h^1)(ka)  ^(J7 ka) 

with  ^(x)  =  xj^(x) ,  fn(x)  2  xh^^(x),  and  e  (o>)  is  the  complex  dielectric 
function  of  the  substrate  sphere.  The  nonradiative  (NR)  rate  can  then  be 
obtained  as  a  difference  between  (1)  and  (2).  On  the  other  hand,  upon 
adopting  the  simpler  nonretarded  IT,  one  obtains  [5] 
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where  the  n-pole  polarizability  ®n(w)  is  given  by 
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Strictly  speaking,  we  cannot  have  radiative  transfer  in  IT  since  the  dipole 
cannot  radiate  in  the  static  limit  [7J.  Nevertheless,  Gersten  and  Nitzan 
(GN)  [2]  have  devised  a  simple  model  to  include  the  radiative  transfer 
within  IT,  and  their  result  for  a  radial  dipole  is  obtained  as 


Distribution/ 
Availability  Codes 
Avail  and/or 


jsot  ]  Special 


L 


□  □ 


(6) 


Zl_ 

70 


|1  +  2( 


2 


Previously,  we  have  studied  in  detail  the  different  distance  behaviors  of 
expressions  (1)  and  (2)  as  compared  to  those  of  (4)  and  (6)  [5-7],  Here  we 
want  to  compare  the  two  theories  by  studying  the  decay  rate  spectrum  as  a 
function  of  the  emission  frequency  w.  Before  we  present  our  numerical 
results,  it  is  interesting  to  point  out  that  the  "resonance  structure" 
implied  from  ET  is  very  different  from  that  obtained  wihtin  IT .  While  in  ET 
resonance  arises  from  the  singular  behavior  of  which  depends  on  both  the 
dielectric  («)  and  geometrical  (a)  properties  of  the  sphere  (the  "Mie 
resonances"),  the  situation  in  IT  is  quite  different  and  depends  only  on  the 
dielectric  property  of  the  sphere.  We  shall  illustrate  these  features 
through  numerical  computations  in  the  next  section. 


3 .  Numerical  Results 

Figure  1  shows  the  total  decay  rates  according  to  both  ET  and  IT  for  a 
radial  dipole  located  at  (d,0,0)  with  d  -  1000  A  above  a  silver  sphere  for 
three  different  values  of  the  radius  a  (with  the  actual  molecule -surface 
distance  d-a  ranging  from  50  to  500  A) .  It  can  be  seen  that  IT  (dotted 
curves)  becomes  more  inaccurate  as  the  molecule  is  located  farther  from  the 
sphere  (note  the  log  scale  for  7j/7q)  [5].  But  what  is  of  more  interest  is 
the  appearance  of  an  extra  resonance  structure  in  ET  (solid  curves)  in 
addition  to  the  ordinary  surface  plasmon  peak  with  ug^  ~  wpU^V/3  [3).  At 
the  low-freqnecy  end,  this  extra  peak  appears  in  ET  with  its  position  and 
magnitude  vary  according  to  the  size  of  the  sphere,  while  only  the  ordinary 
plasmon  peak  shows  up  in  IT.  While  Ruppin  [3J  explained  the  origin  of  these 


peaks  in  terms  of  the  dipole  (I  -  1)  and  higher  multipole  contributions 
(i  >  1),  here  we  present  a  slightly  different  interpretation.  In  fig.  2,  we 
show  the  separate  contributions  from  the  radiative  (dotted)  and  r.onradiative 
(dashed)  transfers  to  the  total  rate  according  to  ET.  It  is  clear  that  this 
new  peak  arises  from  the  resonance  behavior  of  the  radiative  transfer  which 
of  course  is  absent  in  the  nonretarded  IT.  This  peak  may  grow  to  overshadow 
the  ordinary  plasmon  peak,  to  appear  like  a  peak*shift  in  the  decay  rate 
spectrum  (fig.  lc) .  A  similar  peak-shift  phenomenon  has  been  observed 
before  in  surface  enhanced  Rahman  scattering  [8].  In  fig.  3  we  show  the 
same  calculation  for  the  configuration  of  fig.  2  using  IT  and  the  GN  model 
for  the  radiative  rate  (eq.  (6)).  From  the  comparison  of  figs  2  and  3,  we 
see  clearly  that  IT  describes  the  nonradiative  transfers  quite  accurately 

and  the  GN  model  for  the  radiative  rate  within  IT  is  inaccurate. 

IT 

Furthermore,  the  peak  for  7^  according  to  the  GN  model  cannot  move 
according  to  different  sizes  of  the  sphere,  since  it  is  completely 
determined  by  the  dielectric  property  of  the  sphere  as  revealed  in  eq.  (6). 
This,  of  course,  contradicts  the  exact  result  as  shown  in  fig.  1. 

Since  the  exact  dynamics  for  any  dipole-substrate  system  generally 
presents  a  rather  complicated  mathematical  problem,  it  has  been  a  wide 
belief  [2,9,10]  that  under  ordinary  experimental  situations  where  the 
molecule -surface  distance  (s  -  10  A)  is  much  smaller  than  the  emission 
wavelength  ( X  -  10  A) ,  the  simpler  nonretarded  theory  (IT)  can  be  applied 
instead.  Previously,  we  have  pointed  out  that  this  is  only  true  if  the 
conductivity  of  the  substrate  is  relatively  low  so  that  s  is  much  smaller 
than  5,  the  skin  depth  of  the  substrate  [11].  In  this  present  work,  we  h»v* 


observed  further  that  for  a  rough  surface  (a  sphere  in  this  case) ,  the 
condition  s  «  X  is  not  sufficient  if  the  emission  frequency  is  within  the 
resonance  peak  of  the  radiative  transfer  (see  fig.  la  where  s  -  50  A  «  i  - 
10  A) .  The  decay  rate  spectrum  for  a  molecule  near  a  rough  surface 
according  to  the  dynamic  theory  should  have  a  multipeak  structure  in 
general,  and  a  "shift"  in  the  ordinary  peak  may  seem  to  occur  as  explained 
above.  This  should  shed  some  light  on  the  earlier  controversy  as  to  whether 
surface  roughness  may  lead  to  a  peak- shift  in  the  energy  loss  spectrum  as 
observed  in  an  experiment  of  electron  scattering  upon  a  Mo  surface  [12] . 

This  multipeak  phenomenon  does  not  happen  for  the  case  of  a  flat  surface 
(where  the  conditions  s  «  X  and  s  «  6  are  sufficient  for  IT  to  be  valid) , 
since  the  surface  plasmon  mode  cannot  couple  to  the  radiative  mode  without 
the  presence  of  roughness.  However,  the  phenomenon  of  the  morphology- 
dependent  peaks  should  be  very  general  and  should  occur  for  any  kind  of 
rough  surfaces.  In  fact,  we  have  recently  observed  a  similar  multipeak 
phenomenon  for  the  case  of  a  grating  surface  [13] .  Hence  we  conclude  that 
the  exact  dynamic  theory  for  the  classical  decay  rates  of  molecules  at  rough 
surfaces  contains  much  richer  features  which  the  simpler  image  theory  (or 
the  GN  type  approach  to  radiative  decay)  cannot  adequately  describe, 
irrespect  of  the  usual  situation  where  the  molecule -surface  distance  is  much 
smaller  than  the  emission  wavelength. 
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Decay  race  spectrum  for  a  molecule  located  at  d  —  1000  A  from  the  center 
of  a  Ag  sphere  for  various  values  of  the  radius  a.  Curve  a:  a  -  950  A; 
curve  b:  a  -  750  A;  curve  c:  a  -  500  A.  The  solid  curves  represent  the 
results  from  the  dynamic  theory  (ET)  and  the  dotted  curves  the  results 
from  the  image  theory  (IT). 

Decay  rate  spectrum  for  a  molecule  located  at  d  -  700  A  from  the  center 
of  a  Ag  sphere  of  radius  a  -  500  A  according  to  the  dynamic  theorv.  The 
solid,  dotted  and  dashed  curves  denote  the  total,  radiative  and 
nonradiative  rates,  respectively. 

Same  as  Fig.  2,  except  that  the  image  theory  is  applied.  The  solid  and 
dotted  curves  denote  the  total  and  radiative  rates,  respectively. 
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